Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



o 
o 



IRAS 23385+6053: a candidate protostellar massive object * 

F. Fontani 1 , R. Cesaroni 2 , L. Testi 2 , CM. Walmsley 2 , S. Molinari 3 , R. Neri 4 , D. Shepherd 5 , J. Brand 6 , F. 

Palla 2 , and Q. Zhang 7 



o 



o 



Dipartimento di Astronomia e Fisica dello spazio, Largo E. Fermi 2, F50125 Firenze, Italy 

2 INAF, Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 1-50125 Firenze, Italy 

3 IFSI, CNR, Via Fosso del Cavaliere, 1-00133 Roma, Italy 

4 Institut de Radio Astronomie Milimetrique, 300 Rue de la Piscine, F-38406 St. Martin de Heres, France 

5 National Radio Astronomy Observatory, P.O. Box O, Socorro, NM 87801 
VO 6 Istituto di Radioastronomia, CNR, Via Gobetti 101, 40129 Bologna, Italy 

7 Harvard Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

Received date; accepted date 

o . 

i Abstract. We present the results of a multi-line and continuum study towards the source IRAS 23385+6053 

■ performed with the IRAM-30m telescope, the Plateau de Bure Interferometer, the Very Large Array Interferometer 

and the James Clerk Maxwell Telescope. We have obtained single-dish maps in the C ls O (1-0), C 17 (1-0) and 
(2-1) rotational lines, interferometric maps in the CH3C2H (13-12) line, NH3 (1,1) and (2,2) inversion transitions, 
and single-pointing observations of the CH3C2H (6-5), (8-7) and (13-12) rotational lines. The new results confirm 
our earlier findings, namely that IRAS 23385+6053 is a good candidate high-mass protostellar object, precursor 
of an ultracompact Hn region. The source is roughly composed of two regions: a molecular core ~ 0.03 + 0.04 
pc in size, with a temperature of ~ 40 K and an H2 volume density of the order of 10 7 cm -3 , and an extended 
halo of diameter <0.4 pc, with an average kinetic temperature of ~ 15 K and H2 volume density of the order 
of 10 5 cm -3 . The core temperature is much smaller than what is typically found in molecular cores of the same 
diameter surrounding massive ZAMS stars. From the continuum spectrum we deduce that the core luminosity is 
between 150 and 1.6 x 10 4 I/q, and we believe that the upper limit is near the "true" source luminosity. Moreover, 
by comparing the H2 volume density obtained at different radii from the IRAS source, we find that the halo has 
a density profile of the type uh 2 oc r -2 ' 3 . This suggests that the source is gravitationally unstable. The latter 
hypothesis is also supported by a low virial-to-gas mass ratio (Mvm/M gas <0.3). Finally, we demonstrate that 
the temperature at the core surface is consistent with a core luminosity of 10 3 Lq and conclude that we might 
be observing a protostar still accreting material from its parental cloud, whose mass at present is ~ GMq. 

Key words. Stars: formation - Radio lines: ISM - ISM: molecules, continuum - ISM: individual objects: IRAS 
23385+6053 



1. Introduction 

The study of massive stars (M > 10 Mq) and their for- 
mation is important for a better understanding of the 
evolution and morphology of the Galaxy. However, up to 
now most progress has been done in the study of the for- 
mation of low-mass stars (M < 1 Mq). This is a conse- 
quence of the many observational problems which hinder 
the study of the high-mass star formation process: massive 
stars are more distant than low-mass ones, interact more 
strongly with their environment, have shorter evolutionary 
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timescales and mainly form in clusters. Neverthless, re- 
cently a major observational effort has been made to iden- 
tify the very earliest stages of their evolution. Successful 
results have been obtained searching for high-density and 
high-temperature tracers (e.g. NH3 or CH3CN) towards 
selected targets associated with regions of massive star 
formation, such as ultracompact Hn regions (Cesaroni et 
al. Olmi et al. HUM Cesaroni et al. H^Sjl . H 2 

masers (Codella et al. HWI Plume et al. ITW7I Beuther 
et al. 12002 bjl . and IRAS sources (Molinari et al. HMFjI 
Beuther et al. |E002a). Probably the most relevant find- 
ing of these studies is the detection of hot (^100 K), dense 
(<; 10 7 cm -3 ), molecular cores where high-mass stars have 
recently formed. It is believed that such "hot cores" (HCs, 
hereafter) represent the natal environment of O-B stars. 



2 



Fontani et al.: IRAS 23385+6053 



The next step in the search for very young massive 
stars is the detection of a genuine example of massive pro- 
tostar. Molinari et al. £1998b) suggested that the source 
IRAS 23385+6053 might be such an object. 

IRAS 23385+6053 is one of a sample of IRAS sources 
selected in a long-standing project as a likely candi- 
date massive protostar. The selection criteria used in this 
project are based on the IRAS colours (Paha et al. [1991) 
and are aimed at identifying luminous (and hence likely 
massive) stellar objects in the earliest stages of their evolu- 
tion. Observations of molecular tracers such as H2O (Palla 
et al. I1991f> and NH 3 (Molinari et al. 11996(1 have proven 
these objects to be associated with dense molecular gas. 
Additional observations in the continuum with the James 
Clerk Maxwell Telescope (JCMT, Molinari et al. IMMljl 
and the Very Large Array (VLA) (Molinari et al. [1998a) 
have confirmed that a well defined sub-sample of these 
sources is embedded in dense dusty clumps (detected at 
mm and sub-mm wavelengths), and do not present any 
free-free emission, i.e. are undetected at centimeter wave- 
lengths down to levels below that expected on the basis 
of their luminosity. The latter finding suggests that, al- 
though luminous (i.e. massive), the objects embedded in 
the clumps are still too young to develop an Hii region. 
Such a sample contains excellent targets to search for mas- 
sive protostars. 

Recently, Brand et al. ( 200TJ have performed a single- 
dish multi-line study mapping a small number of these 
candidates in various molecular transitions, finding that 
they are associated with clumps that are larger, cooler, 
more massive and less turbulent than those associated 
with ultracompact Hn regions. 

One of these sources, IRAS 23385+6053 (at a kine- 
matic distance of ~ 4.9 kpc), was observed in more de- 
tail at different angular resolutions with the OVRO in- 
terferometer (Molinari et al. 11998b! |2"UT)2"|) and the VLA 
(Molinari et al. l2002[) . This object has also been observed 
with ISOCAM at 7 and 15 /Ltm. The main findings of these 
first studies are the following: 

— the source presents a compact molecular, dusty core 
with diameter ^0.08 pc and mass ^370 M , centered 
at a(J2000) = 23 h 40 m 54!5, o"(J2000) = 61°10'28."10. 
Hereafter, the name 123385 will be used to identify this 
core; 

— the luminosity derived from the IRAS and JCMT con- 
tinuum measurements for a distance of 4.9 kpc is 
~1.6 10 4 L Q ; 

— no continuum emission from the core is detected at 2 
and 6 cm with the VLA, down to ~0.5 mJy beam -1 
(3a); 

— no continuum emission from the core is detected at 7 
and 15 /im with ISOCAM down to ~6 mJy beam -1 
(3a); 

— a faint, small (<0.1 pc) bipolar outflow with axis 
closely aligned with the line of sight has been seen in 
the SiO, HCO+, 13 CO and CS line wings. 



Molinari et al. Ijl998b|l concluded that 123385 is a 
good example of a luminous, young, deeply embedded 
(proto)stellar source. 

Also, Molinari et al. (2002} discovered two extended 
Hii regions with the VLA at 3.6 cm which, however, do 
not coincide with 123385; rather they seem to overlap with 
a cluster of infrared sources detected with ISOCAM that 
surrounds the molecular peak. The infrared emission is 
coincident with a cluster of stars in near-infrared images 
(Molinari, priv. comm.). In this work, we discuss whether 
123385 is a newly-born B0 star just arrived on the ZAMS 
or if it is a massive protostar still in the main accretion 
phase. For this purpose, the question we must answer is 
whether 123385 derives its luminosity primarly from hy- 
drogen burning or from accretion. We believe that un- 
derstanding the nature of 123385 can be improved by ex- 
amining the physical parameters of the core, in particu- 
lar by measuring the temperature of the associated core. 
In fact, in low-mass objects the temperature seems to be 
lower in the earlier evolutionary stages (Myers & Ladd 
1993). By analogy, if high-mass ZAMS stars are found 
in "hot cores" , then high-mass protostars might lie in- 
side colder cores. Furthermore, a detailed picture of the 
source, from small to large spatial scales, should be very 
useful to our purposes. We present here the results ob- 
tained from observations of the C 18 and C 17 (1-0), 
the C 17 (2-1) and the CH 3 C 2 H (6-5), (8-7) and (13- 
12) lines using the IRAM-30m telescope and the Plateau 
de Bure Interferometer (PdBI). We also mapped the NH3 
(1,1) and (2,2) inversion lines with the VLA. In Sect. 2 
we will describe the observations; in Sects. 3 and 4 we will 
present the observational results and derive the physical 
parameters of the source, respectively; in Sect. 5 we will 
discuss these results, and finally we will draw our conclu- 
sions in Sect. 6. 

2. Observations and data reduction 

2.1. IRAM-30m 

The molecular transitions observed with the IRAM-30m 
telescope are listed in Table 1: here, we also give the fre- 
quencies, the half power beam width of the telescope, 
the total bandwidth used and the spectral resolution. 
The main beam brightness temperature, Tmb, and the 
flux density, F v , are related by the expression ^(Jy) = 
4.9T MB (K). 

2.1.1. C 18 and C 17 

C 18 and C 17 data were obtained on September 5, 2000. 
We simultaneously used two 3 mm receivers centered at 
the frequencies of the C 18 (1-0) and C 17 (1-0), and 
two 1.3 mm receiver both centered at the frequency of the 
C 17 (2-1) line to optimize the signal-to-noise ratio. For 
every line, three 3x3 maps with grid spacings 6, 8 and 10" 
and three 5x5 point maps with the same grid spacing were 
made. All maps are centered at q(J2000) = 23 h 40 m 54!5, 
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(5(J2000) = 61°10'28."10: this position corresponds to the 
sub- mm peak detected by Molinari et al. (|2000|> . Map sam- 
pling allows us to cover a region ~ 40" in size. Pointing 
and receiver alignment were regularly checked, and they 
were found to be accurate to within 2". The data were cal- 
ibrated with the "chopper wheel" technique (see Kutner & 
Urich fn)81[l . The integration time was 2 min per point in 
"wobbler-switching" mode, namely a nutating secondary 
reflector with a beam-throw of 240" in azimuth and a 
phase duration of 2 s. 

The system temperature was ~ 250 K at the frequency 
of the C 17 (1-0) and C 18 (1-0) lines, and it was - 1000 
K in the C 17 (2-1) line. 

2.1.2. CH 3 C 2 H 

CH3C2H data were obtained on August 12, 1999. We si- 
multaneously observed the (6-5), (8-7) and (13-12) rota- 
tional transitions in the 3 mm, 2 mm and 1.3 mm bands 
respectively. Only single-pointing observations were made. 
The data have been calibrated using the "chopper wheel" 
technique, and the observations were performed using the 
"wobbler-switching" mode. The system temperature was 
~ 120 K at 3 mm, ~ 260 K at 2 mm and ~ 350 K at 
1.3 mm. 

2.2. JCMT 

An 850//m continuum image was taken on October 16 1998 
with SCUBA at the JCMT (Holland et al. 1155511 towards 
the postion of the core 123385. The standard 64-points 
jiggle map observing mode was used, with a chop throw 
of 2 arcmin in the SE direction. Atmospheric conditions 
were not excellent, with tii^ghz ~ 0.15. Telescope focus 
and pointing were checked using Uranus and the data were 
calibrated following standard recipes as in the SCUBA 
User Manual (SURF). 
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Fig. 1. C 18 and C 17 spectra towards the central posi- 
tion in the maps, which is the peak of the sub- mm emission 
(Molinari et al. 2000). The dashed line indicates the main 
component arising from the source. The dotted line with 
red-shifted velocity shows the second component. In the 
C 17 (1-0), the vertical lines under the spectrum indicate 
the position of the hyperfine components. 



2.3. Plateau de Bure Interferometer 

All transitions observed with the IRAM-30m telescope 
except for CH3C2H (6-5) and (8-7) have been observed 
also with the IRAM 5-elemcnt array at Plateau dc Bure. 
The observations were carried out in November 1998 and 
April 1999 in the B and C configurations of the array. 
The antennae were equipped with 82-116 GHz and 210- 
245 GHz receivers operating simultaneously. These were 
tuned single side-band at 3 mm and double side-band at 
1.3 mm resulting in typical system temperatures of 200 K 
(in the USB) and 1000 K (in the LSB) at 3 mm and 250 K 
(DSB) at 1.3 mm. The facility correlator was centred at 
112.581 GHz in the USB at 3 mm and at 221.921 GHz 
in the LSB at 1.3 mm. The C 18 O(l-0) and C 17 O(l-0) 
lines were covered using two correlator units of 20 MHz 
bandwidth (with 256 channels); two 160 MHz units were 
suitably placed to obtain a continuum measurement at 
3 mm. The C 17 0(2-l) and CH 3 C 2 H(13-12) if =0,1,2 tran- 



sitions were observed with a 40 MHz bandwidth (256 chan- 
nels), while the other K lines were covered with 160 MHz 
bandwidths (64 channels), which were also used to ob- 
tain a continuum measurement at 1.3 mm. The effec- 
tive spectral resolution is about twice the channel spac- 
ing (see "The New Correlator Description of the PdBI", 
http:/ /www.iram.fr). 

Phase and amplitude calibration were obtained by reg- 
ular observations (every 20 min) of nearby point sources 
(2200+420 and 0059+5808). The bandpass calibration 
was carried out on 3C454.3, while the absolute flux den- 
sity scale was derived from MWC349, regulary monitored 
at IRAM. Continuum images were produced after averag- 
ing line-free channels and then subtracted from the line. 
The resulting maps were then cleaned and channel maps 
for the lines were produced. The synthesized beam sizes 
are listed in Tabled 
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Table 1. Transitions observed with the IRAM-30m telescope and the Plateau de Bure interferometer 



Molecular 


Rest freq. 


HPBW (30m) 


Band / Resolution (30m) 


HPBW (PdBI) 


Band / Resolution (PdBI) 


transition 


(GHz) 


(") 


(MHz) 


(") 


(MHz) 


C i8 (1-0) 


109.782 


23 


70 / 0.078 (~ 0.2 km s -1 ) 


2.3x1.94 


16 / 0.156 (~ 0.4 km s -i ) 


C 17 (1-0) 


112.359 


22 


70 / 0.078 (~ 0.2 km s -1 ) 


1.95x1.45 


17 / 0.156 (~ 0.4 km s -1 ) 


C 17 (2-1) 


224.714 


11 


70 / 0.078 (~ 0.1 km s -1 ) 


0.94x0.75 


25 / 0.312 (~ 0.4 km s -1 ) 


CH 3 C 2 H (6-5) 


102.547(*) 


24 


70 / 0.078 (~ 0.2 km s -1 ) 






CH3C2H (8-7) 


136.728(*) 


18 


70 / 0.078 (~ 0.2 km s -1 ) 






CH3C2H (13-12) 


222.166(*) 


11 


140 / 0.078 (~ 0.1 km s -1 ) 


0.94x0.76 


33 / 0.312 (~ 0.4 km s" 1 ) 



(*) rest frequency of the K = transition 
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Fig. 2. From top to bottom: CH 3 C 2 H (6-5), (8-7) and 
(13-12) IRAM-30m spectra towards the peak of the sub- 
mm emission (Molinari et al. 2000). The numbers under 
the spectra indicate the position of the different K com- 
ponents. The Vlsr is relative to the line with K = 0. 

2.4. Very Large Array 

The NRAO 1 Very Large Array (VLA) ammonia ob- 
servations were performed on August 8, 2000. The 
NH 3 (1,1) and NH 3 (2,2) inversion lines at 23.694496 and 
23.722633 GHz, respectively, were simultaneously ob- 
served with spectral resolutions of 1.2 and 4.9 km s -1 
respectively. The array was used in its most compact con- 

1 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, INC. 



figuration (D), offering baselines from 35 m to 1 km. The 
flux density scale was established by observing 3C286 and 
3C147; the uncertainty is expected to be less than 15%. 
Gain calibration was ensured by frequent observations of 
the point source J0102+584, with a measured flux den- 
sity at the time of the observations of 3.4 Jy. The same 
calibrator was also used for bandpass correction. 

The data were edited and calibrated using the 
Astronomical Image Processing System (AIPS) following 
standard procedures. Imaging and deconvolution was per- 
formed using the IMAGR task and naturally weighting 
the visibilities. The resulting synthesised beam FWHM is 
4.0" x 3.7" with position angle —50°, and the noise level 
in each channel is 4 mJy beam -1 and 3 mJy beam -1 for 
the NH 3 (1,1) and NH 3 (2,2) data respectively. 

2.5. Data reduction and fitting procedures 

The reduction of the IRAM-30m telescope and PdBI data 
has been carried out with the GAG-software developed at 
IRAM and the Observatoire de Grenoble. In the next sub- 
sections we briefly outline the fitting procedures adopted 
to analyze the CH 3 C2H and NH 3 spectra. 

2.5.1. CH3C2H fitting procedure 

CH 3 C2H is a symmetric-top molecule with small dipole 
moment (fi = 0.75 Debye). Its rotational levels are de- 
scribed by two quantum numbers: J, associated with the 
total angular momentum, and K, its projection on the 
symmetry axis. Such a structure entails that for each 
J + l — > J radiative transition, J+l lines with K < J can 
be seen (for a detailed description see Townes & Schawlow 
our observations the bandwidth covers all com- 
ponents for the (6-5) and (8-7) transitions, and up to 
the K = 10 component for the (13-12) line. However, we 
detected only lines up to K = 3. 

In order to compute the line parameters, we have per- 
formed Gaussian fits to the observed spectra assuming 
that all the K components of each J + 1 — > J transi- 
tion arise from the same gas. Hence, they have the same 
LSR velocity and line width. Then, in the fit procedure we 
fixed the line separation in each spectrum to the labora- 
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tory value and we assumed the line widths to be identical. 
In Sect. 21 we will use the line intensities derived from 
this procedure to estimate the kinetic temperature of the 
emitting gas. 

2.5.2. NH3 fitting procedure 

The NH3 (1,1) and (2,2) inversion lines show hyperfine 
structure (see i.e. Townes & Schawlow|T975 ). To take into 
account this structure, we fitted the lines using METHOD 
NH 3 (1,1) and METHOD NH 3 (2,2) of the CLASS pro- 
gram. In this case the fit to the NH3 lines is performed 
assuming that all components have equal excitation tem- 
peratures, that the line separation is fixed at the labo- 
ratory value and that the linewidths are identical. This 
method also gives an estimate of the total optical depth 
of the lines using the intensity ratio between the different 
hyperfine components. 

3. Results 

3.1. IRAMSOm observations 

3.1.1. C 18 and C 17 spectra 

In Fig. 1 we show the C 18 and C 17 spectra taken at the 
central position in the maps, corresponding to the peak of 
the sub- mm emission detected by Molinari et al. ((2000). 
Single-pointing spectra of the CH3C2H (6-5), (8-7) and 
(13-12) lines are shown in Fig. 2. The C 18 (1-0) and 
C 17 (1-0) spectra towards the central position clearly in- 
dicate the presence of two velocity components, centered 
at ~ —50.5 km s _1 and at ~ —47.8 km s _1 , respectively 
(see Fig. 1). The component with higher velocity appears 
very strong especially in the C 18 (1-0) spectrum. We 
have fitted the two components with Gaussians. In par- 
ticular, the C 17 (1-0) line has hyperfine structure: thus, 
we have performed Gaussian fits taking into account this 
structure by using METHOD HFS of the CLASS program. 
This procedure also gives an estimate of the line optical 
depth: we will discuss this point in Sect. 13. 1.51 The line pa- 
rameters of the two components at the peak position are 
given in Table[5] Peak velocity (V\ sr ), full width half max- 
imum (Av) and integrated intensity of the line (J Tmb^v) 
are given in Cols. 3, 4 and 5, respectively. In Col. 2 we 
also give the la rms of the spectra. In the CH3C2H ob- 
servations only the lower velocity component is seen. For 
these lines, we have performed Gaussian fits as explained 
111 

Sect. ESI The results are given in Table 

3.1.2. C 18 and C 17 integrated intensity maps 

Maps of the spectral component at —47.8 km s _1 show 
that the C 18 (1-0) (Fig. OJ) and the C 17 (2-1) lines 
peak at (0",— 15") from the central position of 123385, sug- 
gesting that this emission is unrelated to this source. The 
presence of a secondary source to the South was already 
found by Molinari et al. <|1998b|) in the HCO+(1-0) line. 



Only the component centered at —50.5 km s _1 represents 
the emission arising from 123385: hereafter we shall refer 
to this as the "main" component. The secondary com- 
ponent, centered at —47.8 km s _1 , will not be discussed 
further. The integrated intensity maps of the main com- 
ponent are shown in Figs. |3| and 0] We do not show the 
C 17 (1-0) integrated map because the line is faint over 
the whole map (S/N< 3tr in many points). The first clear 
result is that the emission looks quite flat and without an 
obvious peak towards 123385. This is especially evident in 
the C 17 map. Also, in the C 18 (1-0) and C 17 (2-1) 
maps, the half maximum power contour is elongated in 
the N-S direction: this could be due to the presence of the 
secondary component in the spectra. In fact, it is difficult 
to separate the two components, especially in the south- 
ern part of the map where the overlap of the two lines is 
more pronounced. In an attempt to derive the size of the 
emtting region in each line, we have estimated the angular 
diameter considering only the northern part of the maps 
(that with positive offset in S), which is less affected by 
the secondary component: we find an average diameter of 
the emitting region ~ 22" for the C 18 (1-0) line, and 
~ 18" for the C 17 (2-1) line. After correcting for the 
beam size (see Table 0, one finds a deconvolved angular 
diameter of ~ 18" (corresponding to ~ 0.45 pc) for the 
C 18 (1-0) line and of - 15" (~ 0.30 pc) for the C 17 
(2-1) line. 

In conclusion, from the 30-m maps of 123385 in the 
C 18 and C 17 lines, one cannot see an intensity peak 
towards the nominal position of the core, which instead is 
clearly seen by Molinari et al. (1998b) in different molec- 
ular transitions, as well as in our PdBI maps in the same 
lines, as shown in Sect. 13.31 We will further discuss this 
point in Sect. 15.21 

3.1.3. Optical depth of the lines 

The availability of two rotational transitions of the same 
molecular species (namely C 17 (1-0) and (2-1)) allows 
us to derive a map of the kinetic temperature, as we will 
explain in Sect. 14.31 However, this requires knowledge of 
the optical depth of such lines. Since we have also observed 
the C 18 (1-0) line, it is possible to estimate the opacity 
of the C 17 (1-0) line from the ratio of two spectra, as- 
suming equal temperatures and beam filling factors, and 
an abundance ratio of 3.7 for 18 0/ 17 (see Penzias 119811 
Wilson & Rood 1994). We find an intensity ratio typically 
between 3 and 4, with a mean value of ~ 3.5, which implies 
that the C 17 (1-0) line is optically thin. 

For C 17 0, it is possible to derive the optical depth 
from the measure of the relative intensities of the hyperfine 
components. Unfortunately, in our spectra we are not able 
to resolve the hyperfine structure, hence the optical depth 
derived in this way is not reliable. Concerning the C 17 
(2-1) line, we cannot measure the optical depth directly 
from our observations. However, in LTE conditions, one 
can demonstrate that up to temperatures as high as ~ 40 
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Table 2. Line parameters of CO isotopomers at the peak position (IRAM-30m observations). All line parameters have 
been obtained from a two-Gaussian fit, except for the C 17 (1-0) line, for which we have used METHOD HFS of the 
CLASS program to take into account the hyperfine structure (see text). 



Line 


rms 
(K) 


Vl BI 

(km s" 1 ) 


(km s" 1 ) 


J T M B<iv 
(K km s" 1 ) 


C is O (1-0) 


0.065 


-50.0 


2.38 


3.73 






-47.8 


1.20 


2.59 


c 17 o (1-0) 


0.064 


-49.9 


2.08 


1.18 






-47.5 


1.12 


0.41 


C 17 (2-1) 


0.18 


-49.8 


2.93 


5.06 






-47.3 


0.90 


0.69 



Table 3. CH3C2H line parameters (IRAM-30m observations) 



Line 


rms 


VhSR 


A«l/ 2 


j'Tusdv (K km s" 


- 1 ) 




(K) 


(km s _1 ) 


(km s _1 ) 


K = 1 2 


3 


CH3C2H (6-5) 


0.011 


-50.3 


2.48 


0.37 0.29 0.13 


0.07 


CH3C2H (8-7) 


0.024 


-50.5 


2.64 


0.41 0.40 0.22 


0.13 


CH 3 C 2 H(13-12) 


0.032 


-50.5 


2.57 


0.75 0.71 0.41 


0.30 



K, the optical depth of the (2-1) line is < 3 times that 
of the (1-0) line. Since the maximum value of the opacity 
of the (1-0) line derived from the line ratio is r ~ 0.18, 
the (2-1) line must have an optical depth less than ~ 
0.5. Thus, when in Sect. O 

we will derive the kinetic 
temperature from the ratio between the C 17 (2-1) and 
(1-0) lines, we will assume that the lines are optically thin. 

3.2. JCMT observations 

123385 has been mapped in the sub-mm continuum at 
850/xm with SCUBA, with an angular resolution of 15 " 
(Fig.^J. The intensity profile presents an unresolved com- 
pact central peak, and an extended halo which surrounds 
the peak up to ~ 100 ". The figure also shows the inte- 
grated emission of the main component of the C 18 (1-0) 
line (upper panel of Fig.|3J). We can see that there is a rea- 
sonable match between the 850 /jm emitting region and 
that traced by the C 18 (1-0) line, although the corre- 
sponding peaks are not coincident. 

3.3. Plateau de Bure observations 

We present here the results of molecular line and millime- 
ter continuum observations carried out with the Plateau 
dc Bure Interferometer. 

3.3.1. Molecular lines 

All molecular transitions observed with the IRAM-30m 
telescope have been also imaged with the PdBI, with the 
exception of the CH 3 C 2 H (6-5) and (8-7) lines. In Fig.© 
we show the maps obtained integrating the emission under 
the C 18 (1-0), C 17 (2-1) and CH 3 C 2 H (13-12) lines. 
The secondary component is not detected, except perhaps 



in C 18 (1-0), in which the shape of the red wing of the 
line might be affected by the secondary component. This 
is likely due to the fact that the secondary component is 
more extended, thus its emission has been resolved out. 

We do not show the C 17 (1-0) map because this is 
very noisy, and the signal is undetected up to a 3<r level 
of - 36 mJy beam" 1 . The C 18 (1-0), C 17 (2-1) and 
CH3C2H (13-12) emission is detected with a good signal- 
to-noise ratio. We can clearly see the presence of a central 
compact core coincident with 123385. In the C 17 (2-1) 
and CH3C2H (13-12) lines, the source shows a complex 
morphology, while it looks more simple in the C 18 (1-0) 
line. This is consistent with the different angular resolu- 
tion: at 1.3 mm the beam size is ~ 2 times smaller, allow- 
ing us to better resolve the core structure. In particular, it 
can be noted that the C 17 (2-1) emission is marginally 
depressed at the peak position of the continuum and the 
CH3C2H (13-12) line. This might be suggestive of CO de- 
pletion in the inner region of the core (see discussion in 
Sect. 5.2). 

We have estimated the diameters of the emitting re- 
gions from the half maximum power contours shown in 
Fig. if we assume the source to be Gaussian, we can 
derive the angular diameter after beam deconvolution, as 
given in Col. 2 of Table 

In Fig. □spectra of the C 18 (1-0), C 17 (1-0) and 
(2-1) and CH3C 2 H (13-12) lines are shown, that have 
been obtained integrating the emission arising within the 
4cr level inferred from the maps of Fig. El The availabil- 
ity of IRAM-30m spectra allows a comparison between 
interferometric and single-dish data. This comparison is 
also shown in Fig. [7| where the single-dish spectra have 
been resampled with the same resolution in velocity as the 
PdBI spectra, namely ~ 0.4 km s _1 . The superposition of 
the spectra clearly demonstrates that a large fraction of 
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Fig. 3. Top panel: IRAM-30m map of the main compo- 
nent of the C 18 (1-0) line obtained integrating between 
—52.5 and —49.0 km s" 1 . Contour levels go from ~ 1.5 K 
km s _1 (corresponding to ~ 3<t), to the maximum value of 
-4 K km s -1 in steps of 0.5 K km s _1 (~ la). The thick, 
dashed contour indicates the half maximum power con- 
tour. Bottom panel: same as Top for the secondary com- 
ponent of the C 18 (1-0) line, integrated between —49.0 
and —45.0 km s~ x . Contour levels range from 1 K km s _1 
(~ 3cr) to 3.3 K km s _1 in steps of 0.3 K km s _1 (~ la). 



the extended emission is filtered out by the interferome- 
ter. This is especially evident in the C 18 (1-0) line, for 
which the flux measured with the PdBI is ~ 11 times less 
than those obtained with the IRAM-30m telescope. This 
is also true for the C 17 (2-1) line, for which the intensity 
ratio is ~ 10. In the CH3C2H (13-12) transition this effect 
is smaller as the single-dish flux is ~ 3 times larger. As 
noted earlier, the secondary component is resolved out in 
the PdBI spectra. 

Looking at the C 17 (2-1) spectrum of Fig. □ the 
line profile seen by the interferometer seems to peak at 
a slightly blue-shifted velocity with respect to the 30-m 
spectrum. We believe that this is due to the fact that 
the S/N ratio is poor in the PdBI spectrum, hence the 
peak velocity, and the other line parameters, are affected 
by large uncertainties. In fact, the offest observed (~ 0.5 
km s _1 ) lies within the 3a uncertainty on the line peak 
velocity (~ 0.45 km s -1 ) obtained from a Gaussian fit. 




20 10 -10 -20 
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Fig. 4. Same as Fig.[3]for the main component of the C 17 
(2-1) line, integrated between —52.5 and —49.0 km s -1 . 
Contour levels go from ~ IK km s _1 (~ 3a level), to 3.3 
K km s _1 in steps of la. 
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61°09'40" k i i 

23 h 41 rn 00 ! 55 s 50 s 

R.A.(J2000) 

Fig. 5. SCUBA map of IRAS 23385+6053 at 850 fj,m 
(grey scale). Levels range from 0.3 (3c) to 1.9 by 0.2 Jy 
beam -1 . The solid contours represent the C 18 (1-0) map 
(main component) obtained with the IRAM-30m telescope 
(Fig. OH , with the thick line corresponding to the contour 
at half maximum. 

3.3.2. Millimeter continuum 

With the PdBI we detect a compact central core, as can 
be seen from Fig. [HI where we show the 3 mm and 1.3 mm 
continuum maps. The 3 mm map shows a compact core 
with an elongated shape along the E-W direction. This 
shape is better resolved at 1.3 mm, revealing a secondary 
faint peak of emission offset by ~ +2 " in R.A. from the 
main peak. This secondary peak is not detected in any of 
the molecular tracers. The angular diameters are 1.5" and 
1.1" at 3 mm and 1.3 mm, respectively. Measured flux 
densities F v are 0.0124 and 0.148 Jy at 3 and 1.3 mm, 
respectively, obtained by integrating over the solid angle 
down to the 3a level. The flux density of the secondary 
source at 1.3 mm is ~ 5% of the total, namely 0.008 Jy. 
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Fig. 6. Maps obtained with the PdBI. a: C 18 O(l-0) map 
integrated over the velocity range (—49,-52) km s -1 . 
Contour levels range from 0.05 (~ 3(7 ) to 0.14 by 0.02 Jy 
beam -1 . The thick contour indicates the half maximum 
power contour. The spectral rms is ~ 0.014 Jy beam -1 . 
The grey-scale image shows the 3 mm continuum emission, 
ranging from ~ 1.9 (~ 4cr) to 7.4 by 0.8 mJy beam -1 . 
The ellipse at the bottom right represents the synthe- 
sised beam, b: same as a for the C 17 (2-1) line. The 
map has been integrated over (—47,-52) km s . Contour 
levels range from 0.03 to 0.09 by 0.015 Jy beam -1 . The 
grey-scale image represents the 1.3 mm continuum emis- 
sion, whose countour levels range from 6 to 24 by 2.5 mJy 
beam -1 . The spectral rms is ~ 0.012 Jy beam -1 , c: same 
as b for the CH3C2H (13-12) line, integrated under the 
K=0 and 1 lines. Contour levels range from 0.03 to 0.1 by 
0.01 Jy beam -1 . The spectral rms is ~ 0.013 Jy beam -1 

We can compare the continuum flux density measured 
with the PdBI and with SCUBA, to estimate the amount 




-50 -40 -30 -20 

V lsr (km s ) 

Fig. 7. Flux density comparison between IRAM single- 
dish (dashed line) spectra and Plateau-de-Bure interfero- 
metric spectra (solid line). The interferometric flux den- 
sities have been multiplied by a factor of 10 for the CO 
isotopomers, and by a factor of 3 for the CH3C2H lines. 
The velocity of the CH3C2H spectrum is computed with 
respect to the frequency of the line with K = 0. 



of the flux density filtered out by the interferometer. Using 
a spectral index for the dust emissivity of 4 (see Sect. l4.l|l . 
we have extrapolated the flux density observed at 850 /jm 
with SCUBA, to 1.3 mm. This has been smoothed to an 
angular resolution of ~ 22", corresponding to the primary 
beam of the PdBI at 1.3 mm. The flux density measured 
in this map towards the core 123385 is ~ 0.28 Jy. We 
estimate in ~ 0.11 Jy the flux density to arise from the 
halo: this is filtered out by the interferometer, and thus the 
flux density that we should observe with the PdBI is ~ 
0.17 Jy. Taking into account the uncertainties that affect 
our measurements (in particular the combined calibration 
errors of both SCUBA and PdBI, which are at least of 
the order of 30%), the 0.15 Jy observed at 1.3 mm with 
PdBI is consistent with what is expected. Thus, within 
the errors, it is plausible that with the PdBI we observe 
all the continuum flux arising from the core 123385. 
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Fig. 8. Top panel: map of the 3 mm continuum emission. 
The rms level is ~ 0.6 mJy beam -1 and the signal-to- 
noise ratio is ~ 15. Contour levels range from 1.6 to 7.8 
by 0.6 mJy beam -1 . The thick, dashed contour represents 
the half maximum power contour. Bottom panel:same as 
top for 1 mm continuum emission. The rms is ~ 2 mJy 
beam -1 at 1.3 mm and the signal-to-noise ratio is ~14. 
Contour levels range from 6 to 24 by 2 mJy beam -1 . The 
ellipses in the bottom right are the synthesized beams. 



3.4. Very Large Array A/f/3 maps 

With the VLA, we have mapped the NH 3 (1,1) and (2,2) 
inversion transitions. The corresponding maps and inte- 
grated spectra are shown in Fig. [!|] In both maps, the 
NH3 emission is superimposed on the C 18 (1-0) map. 
As for the CO isotopomers observed with the IRAM-30m 
telescope, the emission does not peak towards the position 
of the core detected at 3 mm and 1.3 mm with the PdBI. 
Although the signal-to-noise ratio is poor in both maps 
(see Fig. , one can see that the emission arises from a 
"clumpy" structure, and it peaks at ~ 5" (~ 0.12 pc) from 
the core position. In Sect. l5.2l we will discuss this result by 
comparing the NH3 maps with those obtained in the other 
tracers. Comparison between Effelsberg-lOOm (Molinari et 
al. 11996ft and VLA observations are shown in Fig. ^] For 
the main component (centered at —50.5 km s -1 ), with the 
VLA we lose about the 50% of the total flux density seen 
with the single-dish telescope; for the NH 3 (1,1) line, this 
occurs for all the hyperfine lines. Finally, we completely 
resolve out the secondary component (centered at —47.8 
km s -1 ) as in the other molecular tracers observed with 
the PdBI . 
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Fig. 10. Flux density comparison between Effelsberg- 
100m (dashed line) and VLA (solid line) spectra of the 
NH3(1,1) and (2,2) lines. VLA fluxes have been multi- 
plied by a factor of 2. The thick vertical lines under the 
NH3(1,1) spectrum indicate the position of the two com- 
ponents. The Effelsberg spectra have been resampled to 
the channel spacing of the VLA spectra. 



4. Derivation of physical parameters 

Our goal is to investigate the hypothesis of Molinari et 
al. (1998b), namely that 123385 is a high-mass protostel- 
lar object, by examining the physical parameters of the 
core. First, we will analyze the continuum spectral en- 
ergy distribution (SED). Then, from the PdBI map of the 
CH3C2H (13-12) line, we will compute the core temper- 
ature; finally, from continuum and line emission we will 
obtain the mass and the H2 volume density of the core. 

4.1. Continuum SED and source luminosity 

In discussing the continuum SED, it is important to dis- 
tinguish between two "sub-regions" : 

- the core 123385; 

— a stellar and diffuse emission from a cluster surround- 
ing the core (see Fig. mi. 

The cluster is detected as a patchy "ring" at FIR wave- 
lengths with ISOCAM (Fig. IT2*|l and it seems to overlap 
with two extended sources detected at 3.6 cm with the 
VLA. From the radio flux Molinari et al. (2002) have es- 
timated the spectral type of the ionizing stars to be B0. 5. 

In Fig.lTalthe continuum SED of 12338 5 is shown. With 
respect to Fig. 2 of Molinari et al. ( 1998b), Fig.ll3lcontains 
more data, namely PdBI, IRAS, JCMT, OVRO, ISOCAM 
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Fig. 11. Image: IRAS 23385+6053 observed at 2 fim with 
the IRCAM of the Palomar telescope in the Ks band 
(Molinari priv. comm.), which shows the cluster of stars 
around the core 123385 (indicated by the cross at map 
center). The dashed contours represent the emission of 
the two extended radio sources detected at 3.6 cm with 
the VLA: contour levels range between 0.3 and 3 by 0.3 
mJy beam" 1 (Molinari et al. l2002|l . 




61°1O'O0" - 



23 h 41 m 00 s 23 h 40 m 55 5 23 h 40 m 50 s 

R.A.(J2000) 

Fig. 12. Image of IRAS 23385 +6053 at 15 /im obtained 
with ISOCAM (Molinari et al. H998bjl . The thick polygon 
represents the 3cr level of the HCO + (1-0) line, and it 
identifies the core region. The dashed polygons indicate 
the two clusters of stars which lie around 123385, called 
"ring" in the text. 



and MSX 2 . Since no 7 and 15 /im emission is detected 
with ISOCAM towards the source, we take the 3<r level as 

2 MSX images have been taken from the on-line MSX 
database http:/ /www. ipac.caltech.edu/ipac/msx/msx. html 




1 2 3 4 5 



Log 1Q (X[ / am]) 

Fig. 13. Continuum SED of 123385. Open hexagons and 
stars indicate MSX and ISOCAM flux densities respec- 
tively measured integrating over the "ring" region of 
Fig. El The filled triangle represents the SCUBA mea- 
surement. Filled squares, open squares, filled circles and 
open circles indicate respectively IRAS, JCMT, PdBI and 
OVRO data (see also Fig. 2 of Molinari et al. H998bjl . The 
arrows on the bottom left and bottom right respectively 
indicate the ISOCAM and VLA upper limits (3cr level) of 
the core flux density. The solid line represents a grey-body 
fit, with dust temperature T = 40 K and dust absorption 
coefficient (3 oc v 19 , while the dotted line is a grey-body 
fit with the same (3 and T — 15 K. The dashed line cor- 
responds to a power-law extrapolation of the flux density 
arising from the "ring" region. 



an upper limit at these wavelengths. The same was done 
for the VLA measurements at centimeter wavelengths. 

We performed a grey-body fit to the continuum SED 
of Fig. El (solid line). The resulting bolometric luminos- 
ity is the same as in Molinari et al. jT998b), namely 
L ~ 1.6 10 4 Lq. The slope of the millimetric part of 
the spectrum is fitted assuming a dust opacity k v — 
fe 3 oGHz(^/230(GHz))' 3 , where fe 30G Hz = 0.005 cm 2 g-\ 
which assumes a gas-to-dust mass ratio of 100 (see 
Preibisch et al. 1993). The best fit is obtained for F v oc 
v 3 - 9 . Therefore, the best fit opacity index is (3 ~ 1.9. The 
other best fit parameters are: source temperature of 40 K, 
diameter of 8"and a mass of ~ 400Af Q . Since most of the 
luminosity is due to the IRAS 60 and 100 flux den- 
sities, it is important to assertain the origin of the IRAS 
flux densities. In fact, the IRAS beam is much larger than 
the source diameter (~ 2 arcmin at 100 ^m), thus the 
IRAS flux densities may arise both from 123385 and from 
the cluster of stars surrounding it (the "ring" in the 15 
/im ISOCAM map of Fig. lT2*|l . Hence, the source luminos- 
ity estimated before (L ~ 1.6 10 4 L©) is to be regarded 
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as an upper limit for the luminosity of 123385. However, 
we know from the JCMT and interferometric maps that 
the flux densities at 850 /mi and at longer wavelenghts cer- 
tainly come from 123385. Thus, we performed a grey-body 
fit using only these data (dotted line of Fig. I13|l : the cor- 
responding grey-body fit gives a luminosity of ~ 150 Lq 
which represents a lower limit. The two limits thus ob- 
tained are very different, but we believe that the source 
luminosity is closer to the upper limit. In fact, the MSX 
image at 21 /jm shows that the emission at this wavelength 
arises from the "ring" , hence very likely this region is also 
responsible for the 25 /im emission. However, as already 
discussed by Molinari et al. (|1998bli . it is unlikely that the 
"ring" contributes significantly to the IRAS flux densities 
at 60 and 100 /im, because even a power-law extrapola- 
tion of the spectrum of the extended region (dashed line 
in Fig. I13|l gives flux densities at 60 and 100 /im much 
lower (a factor of ~ 10) than the values observed. We 
conclude that most of the emission at 60 and 100 /im is 
due to 123385, and therefore a luminosity of ~ 10 4 L 
is not very far from the real value. We will further dis- 
cuss this point in Sect. 15.31 on the basis of our CH3C2H 
observations. 

It is also worth noting that the 1.3 mm flux density 
measured with the PdBI is ~ 3 times less than the value 
measured with the JCMT at the same wavelength. In fact, 
as explained in Sect. 13.3121 with the interferometer we ob- 
serve approximately all the flux density arising from the 
compact core, but we miss that coming from the extended 
halo, which we do detect with the JCMT. Finally, it is un- 
likely that the source associated with the secondary peak 
detected in the 1.3 mm map contributes significantly to 
the total observed luminosity, because its 1.3 mm flux den- 
sity is only 5% of that from the main core. 

4.2. Kinetic temperature and column density from 
CH3C2H rotation diagrams 

From the CH 3 C 2 H (13-12) spectra of Figs. 2 and we 
derive the kinetic temperature and the total column den- 
sity of the molecule by means of the Boltzmann diagrams 
method (see e.g. Fontani et al. 12002)1 . The fundamental as- 
sumption of the method is the local thermodynamic equi- 
librium (LTE) condition of the gas. Such an assumption is 
believed to work very well for CH3C2H because of its low 
dipole moment (see Bergin et al. 1994). Under the further 
assumption of optically thin emission, the column density 
Ni of the upper level i of each transition is proportional 
to the integrated intensity of the line. 

In Fig. ll4I Boltzmann diagrams are shown: the straight 
lines represent least square fits to the data, and the ki- 
netic temperature and the column density are related 
to the slope and the intercepta, respectively. The inter- 
ferometric data yield T k ~ 42 K and AT tot ~ 4 10 14 
cm~ 2 for the core. It is evident that such temperature 
is well below that of a HC. In a recent paper, Thompson 
& Macdonald ( 2003 ) have observed rotational transitions 



of the CH3OH molecule towards IRAS 23385+6053, from 
which they derive that the kinetic temperature is likely 
to be smaller than 50 K, and therefore it is unlikely 
that IRAS 23385+6053 might contain a HC. Our study 
strongly supports this idea and give a much more accu- 
rate temperature estimate of the core. 

With the IRAM-30m data we obtain T k ~ 34 K and 
JVtot — 5 10 15 cm -2 . In this case, the column densities Ni 
have been corrected for the beam filling factor. For this 
purpose, an estimate of the CH3C2H emitting region is 
needed. Unfortunately, we have no direct estimate because 
no CH3C2H single-dish maps are available. However, we 
can use the source size estimated from the SCUBA map, 
given in Sect. 4.5. 

The previous results are summarised in Table ^ The 
errors on Tk and N to t have been computed in two dif- 
ferent ways: for the single-dish data we must take into 
account the calibration error. This error affects simulta- 
neously each line in the same bandwidth. Hence, in the di- 
agrams, the effect of this error is to shift the points in the 
same bandwidth (e.g. the (6-5) lines) by the same quan- 
tity. This error is ~ 10% for the (6-5) lines and ~ 20% 
for the (8-7) and (13-12). Hence, in order to compute the 
uncertainties on and A^t, we have varied the values of 
Ni by 10% simultaneously for all the (6-5) K components, 
and by 20% for the (8-7) and (13-12). These variations 
have been made in different directions. The uncertainties 
in Tableware the maximum difference between the values 
thus obtained and the nominal ones. For the PdBI esti- 
mates, we have data of a single band, and the calibration 
error does not affect Tk but only N to t- 

Let us now discuss the two basic assumptions of the 
method, i.e. optically thin emission and LTE conditions 
of the gas. The LTE assumption holds very well when 
the H2 density is greater than a critical value (see e.g. 
Spitzcr 1978). For the observed transitions the maximum 
critical density is ~ 10 5 cm -3 , and we shall see in Sect. l4.5l 
that the inferred H2 density is above this value, justifying 
the assumption of LTE conditions. Concerning the optical 
depth of the lines, the excellent agreement between the 
data and the linear fit of Fig. [21 strongly supports the 
assumption of optically thin lines. In fact, large optical 
depths are expected to affect mainly the transitions with 
lower excitation. Hence we should have a flattening of the 
Boltzmann diagram at lower energies, which is not seen in 
Fig.Hl 

4.3. Physical parameters from CO isotopomers 

In the previous section we have computed the kinetic tem- 
perature of the source from the CH3C2H lines. We can also 
estimate the temperature by means of the C 17 lines ob- 
served with the IRAM-30m telescope. Using a source size 
of 15", found in Sect. EH~21 from the C 17 (2-1) map, we 
have derived the beam-averaged brightness temperatures 
for the C 17 (1-0) and (2-1) lines. Then, from the inte- 
grated line intensities, and assuming optically thin lines 
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and LTE conditions, wc find a "Boltzmann relationship" 
between the total column density of the upper level of each 
transition and the corresponding rotational energy, as dis- 
cussed above for the CH3C2H lines. We have done this for 
each point of the C 17 maps in which the C 17 (1-0) line 
is stronger than the 3<r level. The resulting average ki- 
netic temperature is ~ 15 K. It is important to stress that 
in deriving the kinetic temperature from C 17 line ratios 
we have assumed LTE conditions. Statistical equilibrium 
calculations (see Wvrowski ll997fl show that this assump- 
tion holds if «h 2 > 10 5 cm~ 3 . Otherwise, the temperature 
estimated in this way is to be regarded as a lower limit. 
Since we obtain njj 2 ~ 10 5 cm~ 3 over the whole map, we 
conclude that the LTE assumption is satisfied in our case. 

We obtain a total column density ranging from 9 10 
to 7 10 15 cm -2 : the C 17 total column density given in 
Table 0] namely 4 10 15 cm -2 , is the mean value over the 
map. 

4.4. Physical parameters from NH 3 (1,1) and (2,2) 

The NH3 emission seen with the VLA arises from a re- 
gion whose diameter is "intermediate" between that of the 
core 123385 and that of the C 17 and C 18 total emitting 
region seen with the IRAM-30m telescope. The emission 
peak does not coincide with the continuum nor with the 
C 18 (1-0) line. The NH 3 (1,1) hyperfine structure allows 
us to estimate the optical depth as explained in Sect. 12.51 
we obtain r <C 1 for this line. For the (2,2), we have no 
direct estimates, but we can reasonably assume that for 
temperatures between ~ 10 and ~ 40 K, under LTE condi- 
tions, the (2,2) line is also optically thin. The temperature 
and column density can be derived following the method 
by Ungerechts et al. (|1986(l . We obtain a temperature of 
~ 26 K and a total column density of 6 10 15 cm -2 . We 
have assumed equal linewidth and beam filling factor for 
both transitions. The results are given in Table 01 

4.5. Mass and density estimate 

The grey-body fit to the spectrum of Fig. El is obtained 
for a dust temperature of 40 K, a dust opacity k v oc v ls> , 
an average source diameter of 8", (~ 0.2 pc), and a mass 
of - 400 M Q . 

It is also useful to derive the mass of the core seen 
with the PdBI. From the 3 mm continuum flux density, 
using the core temperature deduced from the CH3C2H 
lines (namely T = 42 K), fc„ cx v 19 , and assuming opti- 
cally thin dust emission, we find that the mass contained 
inside ~1."5 is ~ 150 M Q . Note that we have used the 
3 mm flux density, rather than that at 1.3 mm, because 
we expect to miss less flux at this wavelength. 

To compute the mass from the molecular line emission, 
we use the expression: 

7rZ)(cm) 2 iVt t(cm- 2 )mH 2 (g) 



where D is the source diameter, toh 2 the mass of the 
H2 molecule, X the abundance of the molecule relative 
to H2, and -/Vtot is the source averaged column density. 
X = 1.2 10" 7 , X = 3.4 10~ 8 and X = 2 10~ 9 for the 
C 18 0, C 17 and CH3C2H mean abundances, respectively. 
The relative abundances of the CO isotopomers have been 
computed from Wilson & Rood l|1994f) for a galactocentric 
source distance of 11 kpc. The CH3C2H abundance is an 
average value from Fontani et al. (2002) and Wang et al. 
( 1993), who studied the CH3C2H emission in massive star 
formation sites. The resulting values of Mn nc are reported 
in Col. 7 of Table 5. 

From the line-widths of the observed transitions and 
the corresponding angular diameters in Table we can 
also derive the mass required for virial equilibrium: assum- 
ing the source to be spherical and homogeneous, neglect- 
ing contributions from magnetic field and surface pressure, 
the virial mass is given by (MacLaren et al.[1988): 

M V ir(M ) = 0.509 d(kpc) 9 s (arcsec) Au 2 /2 (km/s) (2) 

In Table[S]we show the resulting virial mass deduced from 
the CH3C2H (13-12) lines observed with the PdBI and 
also from other lines. 

Finally, from Mi; ne we derive the H2 volume density, 
fiH 2 . All masses and densities are given in Tables |S] and [fjj 
It is worth noting that the density of the central core is 
> 10 7 cm~ 3 : this value is well above the maximum critical 
density of the CH 3 C2H lines observed, thus supporting our 
LTE assumption and our derivation of the kinetic temper- 
ature. As usual, the main errors in estimating the masses 
are due to the molecular abundances, the assumed dust 
opacity and the gas-to-dust ratio, and it is very difficult 
to quantify the uncertainties for these parameters. 

One can see clearly that masses and densities derived 
from different tracers differ very much from one another. 
This is partly due to the fact that different tracers arise 
from different regions. However, even for the same re- 
gion, the mass derived from the continuum emission is 
a few times that estimated from the line emission. This 
is likely due to uncertainties in the assumed molecular 
abundances. We shall come back to this point in Sect. 15.21 

5. Discussion 

The availability of high and low angular resolution obser- 
vations of IRAS 23385+6053 allows us to paint a detailed 
picture of the source, and discuss the nature of the core 
123385. 

5.1. Virial equilibrium and density structure 

The values of the masses estimated for the same region 
using dust and line emission show a systematic difference: 
from Table and |H1 one can see that the mass estimated 
from dust (M cont ) is ~ 2-3 times greater than that deduced 
from molecular lines (Mi; ne ). However, Mii ne is affected by 
several problems which may lead one to underestimate the 
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Table 4. Temperature and column density estimates. 9 S and D are the angular and linear source diameter, 
respectively.Tk is the kinetic temperature, and iV to t the total column density of the corresponding molecular species 



Tracer 


0s 


D 


T k 


iVtot 




(") 


(pc) 


(K) 


(cnr 2 ) 


CH3C2H (PdBI) 


1.4 


0.033 


42±1 


(1.8 ±0.2) 10 14 


NH 3 (VLA) 


~ 5 


0.12 


26±6 


~ 6 ± 2 10 15 


CH3C2H (30m)* 


18 


0.43 


34±6 


(5 ±0.2) 10 15 


C 17 (30m)* 


15 


0.35 


15±5 


~ 4±2 10 15 



* corrected for source size 

Table 5. Mass and H.% column- and volume- density estimates from lines. Molecular abundances relative to H 2 have 
been assumed equal to 1.2 1(T 7 , 3.4 1(T 8 and 2 1CT 9 for C 18 0, C 17 and CH 3 C 2 H, respectively. Nr 2 and Av 1/2 are 
the H 2 column density and the FWHM of the lines, respectively 



Tracer 




D 


Nh 2 




Mvir 


M linc 


nn 2 




(") 


(PC) 


(cm" 2 ) 


(km s" 1 ) 


(M e ) 


(M Q ) 


(cm" 3 ) 


C ly O (2-1) (PdBI) 


1.3 


0.03 


5.4 10 24 


2.9 


30 


42 


6.0 10' 


CH3C2H (13-12) (PdBI) 


1.4 


0.033 


3.3 10 24 


3.0 


32 


15 


1.7 10 7 


C 18 (1-0) (PdBI) 


1.9 


0.045 


4.3 10 24 


3.4 


55 


69 


3.2 10 7 


HCO+(l-0)( t ) 


7 


0.16 


1.66 10 24 


3.3 


200 


370 


3.4 10 6 


C 17 (2-1) (30m) 


15 


0.35 


1.1 10 23 


2.9 


314 


110 


1.0 10 5 


C 18 (1-0) (30m) 


18 


0.43 


1.0 10 23 


2.4 


259 


200 


1.1 10 5 



( f ) From Molinari et al. ( 2002) 

Table 6. Mass and H 2 column- and volume- density estimated from continuum observations 



Tracer 


e s 


D 


Nh 3 


M cont 


n H2 




(") 


(PC) 


(cm" 2 ) 


(Me) 


(cm" 3 ) 


3 mm (PdBI) 


1.5 


0.036 


1.5 10 2b 


150 


14.0 10 Y 


3.4 mm (OVRO f ) 


4 


0.09 


0.4 10 25 


324 


1.6 10 7 


SED (*) 


8 


0.20 


1.5 10 24 


400 


2.5 10 6 


850 Aim (SCUBA) 


18 


0.43 


3.6 10 23 


520 


2.8 10 5 



t From Molinari et al. £l998b I 
* From the SED of Fig. [HI 



mass, as we will discuss in Sect. 15.21 On the other hand, 
the dust emission is independent of molecular abundances, 
so that the mass estimated from the mm continuum seems 
to be more reliable. Since the source is a candidate massive 
protostar, it is useful to discuss its virial equilibrium. In 
order to do this, we need an estimate of the virial mass. We 
have used the diameters deduced from each tracer to give 
an estimate of the virial mass. Then, we have compared 
these estimates to those deduced from the column density 
derived from the same tracer, and with that obtained from 
the continuum emission arising from the same region. The 
result is shown in Fig. The empty circles indicate the 
estimates derived from lines, while the full circles represent 
the continuum. 

In a cloud with a density profile which follows a power 
law of the type n(r) oc r~ m , one can demonstrate that 
the mass contained inside a given radius r, M(r), is 
proportional to r 3 ~ m . Instead, the virial mass depends 
on the source diameter and on the linewidth. However, 
based on the values in Table we conclude that the 



linewidth does not change significantly at different radii. 
Therefore, the estimate of the ratio M(r)/MviR scales ap- 
proximately as r 2 - m . For IRAS 23385±6053, as we shall 
show later (see Fig. llfiJ) . we derive that m=2.3±0.2, hence 
M/Mvm oc r-°- 3 . In Fig. EH we plot the ratio M/M V m 
as a function of D. The curve is in rough agreement with 
data points obtained from the continuum measurements. 
However, the "true" virial mass, deduced from the virial 
theorem, must be computed from the radius of the whole 
cloud and the corresponding Avi/ 2 . We do not know the 
"total" dimension of the source, but we can compute at 
which radius the ratio M(r) /Mvir becomes equal to 1. For 
a power-law density distribution of the type n(r) oc r~~ m , 
the virial mass obtained from Eq. J2J) must be multiplied 
by a factor 1 5 3 ~ 2 " (see MacLaren et al. 119880 , which is 
~ 0.35 for m=2.3. This increases the ratio M/Mvir by a 
factor 3, implying that M/Mym > 3 at a radius of 10", 
and that the virial mass equals the gas mass at a radius of 
~ 50 ", which is much greater than the maximum diameter 
observed by us. Hence, since the real source size is likely 
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Fig. 14. Top panel: rotation diagram inferred from 
CH3C2H lines observed with IRAM-30m telescope. Filled 
circles, filled triangles and open circles indicate respec- 
tively the (6-5), (8-7) and (13-12) transitions. The 
straight line represents a least square fit to the data. 
Column densities are source-averaged values using a 
source diameter of 15" (see Sect. l4.2)l . Bottom panel: same 
as Top for CH 3 C 2 H (13-12) observed with PdBI. The as- 
sumed source size is 1.4". The open circle represents the 
3cr level of the spectrum, taken as an upper limit for the 
K = 4 component. 



to be smaller than that corresponding to M/Mvm = 1, 
we believe that 123385 is likely to be gravitationally un- 
stable. This conclusion is further supported by discussing 
the density profile. 

In Fig.^|we plot the density estimated from all tracers 
observed against the corresponding linear diameter (den- 
sities and diameters are given in Table [3] and |BJ). The data 
follow the relation nn 2 oc D -23 , hence nn 2 oc r~ 23 . This 
is similar to that predicted by star formation models (Shu 
et al. I1987[l , in which star formation occurs in molecular 
clouds which are singular isothermal spheres with density 
oc r~ 2 . This sphere rapidly undergoes inside-out collapse. 
Then, the density profile in Fig. 1161 supports the idea that 
IRAS 23385+6053 is unstable, as expected for a protostar 
which is still in the accretion phase. 



0= M line /M vi: 
•= M cont /M v 



C) 



o 



o 



I 



D (pc) 

Fig. 15. Ratio between masses deduced from gas (My lnc ) 
or dust (Af cont ) and the corresponding virial mass as a 
function of the angular diameter. All values have been 
taken from Tables and H3 Filled circles assume masses 
of the region estimated from dust, while empty circles as- 
sume gas masses from the total column density of observed 
lines. The dotted line represents the ratio M/Mvm> ex- 
pected to be oc r~ - 3 across the source. The estimated 
errors are of the order of ~ 30%. 




850 (ira (SCUBA) 



Log 10 [D (pc)] 

Fig. 16. H2 volume density against the linear diameter D 
for all the tracers observed. All values have been taken 
from Table The straight line is a linear fit to the data, 
which gives nn 2 oc D -2 ' 3 . 

5.2. A depleted core? 

On the basis of Tables 4-6 and the SCUBA map of Fig-EJ 
one can identify two regions in 123385: a compact core and 
a surrounding extended halo. Their physical parameters 
are: 

- Di = 0.035 pc, Ti ~ 42 K, m ~ 14.0 10 7 cm" 3 for the 
core; 

- D 2 ~ 0.4 pc, T 2 ~ 15 K, n 2 ~ 3 10 5 cm~ 3 for the halo. 

where the diameters and njj 2 volume densities have been 
derived from the continuum measurements, and the ki- 
netic temperatures from the lines. Now, using a simple 
source model, we want to understand if it is possible to 
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clarify the absence of a central peak in single-dish maps, 
and also the non- detection of C 17 (1-0) in the PdBI ob- 
servations. 

Let us consider a source consisting of the two regions 
outlined above: a core and a surrounding halo. Both are as- 
sumed to be spherical, homogeneous and isothermal. The 
expected brightness temperature Tb of the C 17 (1-0) 
line for such a source is shown in Fig. 1171 where Tb is 
plotted against the radial distance R from the source cen- 
ter. Convolving Tb with a gaussian beam, we derive the 
corresponding main beam brightness temperature Tmb- 
Fig.^|represents Tb convolved with a gaussian beam with 
HPBW=22" (that of the IRAM 30-m telescope); one can 
see that the emission coming from the core is flattened by 
beam dilution effects. 

We now discuss why the central core is not detected in 
the C 17 (1-0) PdBI map. Since the molecular emission 
is extended, we must determine first how much flux den- 
sity is filtered out by the interferometer. Starting from the 
brightness temperature profile of Fig. 1171 we have used a 
modified version of the procedure described in Wilner et 
al. 1|2000[) and Testi et al. <|2DU3Jl . The program uses the 
model to create first the source image in the sky; then, it 
computes the Fourier Transform and resamples the visibil- 
ity function in the UV points corresponding to those mea- 
sured in that channel, which finally are used to reproduce 
the image. In order to compare the PdBI map with what 
is predicted by the model, we plot in Fig. ^| the channel 
where the peak is observed (at —50.5 km s" 1 ), superim- 
posed on the map predicted by our model. We have also 
done the same using a model that reproduces the "halo" 
only. It is evident that the model reproduces a central peak 
which is instead completely lost in the PdBI map; on the 
other hand, the emission of the halo is quite consistent 
with what observed if we consider that our model is very 
simple, and does not consider possible patchy structures 
inside the halo. 

A possible explanation of the discrepancy between the 
PdBI map and our model can be that in the core the C 18 
and C 17 species are partially depleted because they are 
frozen onto dust grains. This is suggested by the fact that 
gas masses deduced from the CO isotopomers given in 
Table are systematically lower than the dust masses 
inferred from the continuum measurements. This might 
be due to a partial freeze-out of the molecules on the 
dust grains. This situation is possible even at tempera- 
tures as high as 40 K if dust grains have ice mantles. In 
fact, chemical models predict that ice mantles may exist 
at the densities and temperatures measured by us in the 
core (Van Dishoeck & Blake [19983, and this allows par- 
tial depletion of CO and its isotopomers. Furthermore, 
recently Thompson & Macdonald ( 2003 ) have performed 
a chemical analysis of 123385 using molecular lines with 
rest frequencies in the range 330-360 GHz, finding that the 
chemical composition of 123385 is consistent with that of a 
molecular core in the middle evaporation phase, i.e. when 
the majority of the molecular species are beginning to be 
evaporated from dust grains ice mantles. 



C O (1-0) 



4 6 
R (arcsec) 



Fig. 17. Brightness temperature of the C 17 (1-0) line 
computed assuming a source as described in the text, and 
a C 17 relative abundance of 1.4 10~ 8 . R represents the 
radial distance from the source center. 



c 17 o (1-0) 

HPBM = 22 orcsec 



2 4 6 8 

R (arcsec) 

Fig. 18. Brightness temperature of FigHTI convolved with 
the IRAM-30m telescope beam (22"). 



Also, dust grains with ice mantles might be responsible 
for the non-detection of NH3 towards the source center: 
in fact, the same chemical models also predict that the 
presence of ice mantles allows NH3 depletion up to tem- 
peratures of 90 K. Finally, it is worth noting that partial 
depletion of CO is consistent with detection of CH3C2H 
emission in the core. In fact, Ruffle et al. (1997) found 
that in dense cores, when CO depletion occurs, molecu- 
lar species produced from CH and without oxygen (like 
CH3C2H) increase their rate of production. 



5.3. The nature of 123385 

Here we want to find out whether our results may be used 
to establish the nature and the evolutionary stage of the 
embedded YSO. 

In cores where the gas is heated by an embedded star, 
the gas temperature at a radial distance R — 6d from the 
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central star is expected to scale as (see e.g. Plume et al. 

Tcx (^r) < 3 > 

where L is the luminosity of the embedded star. Eq. (|3J 
represents the gas temperature, but at densities as high 
as 10 6 cm~ 3 we can assume coupling between gas and 
dust. Hence, at such densities Eq. © also describes the 
dust temperature. The index a depends on the opacity 
index, /?, (Doty & Leung ITT)94[1 . and it lies between 0.3 
and 0.5. The steepest power law is obtained when a = 
0.5, and it corresponds to the optically thick case, when 
the luminosity irradiated at the core surface R equals the 
stellar luminosity, namely: 

L* = 4irR 2 aT 4 , (4) 

where T represents the surface temperature. If instead the 
core is optically thin, the temperature profile is flatter. 
Therefore, the core temperature predicted by Eq. rep- 
resents a "lower limit", and at 0.017 pc (corresponding 
to the radius of 123385 deduced from the CH 3 C 2 H PdBI 
emission) for a central star with luminosity L ~ 1.6 10 Lq 
this lower limit is ~ 51 K. The temperature derived from 
CH3C2H is significantly lower than that. This discrepancy 
is confirmed by theoretical models: Osorio & Lizano ( 1999 ) 
have performed a theoretical simulation of the temper- 
ature profile of 123385 assuming a source luminosity of 
~ 1O 4 L0, deriving a temperature much higher than 42 K 
at the same radius (see their Fig. 5). To compare our result 
with what found in HCs, in Fig. [2U| we plot the temper- 
ature T at the core surface versus L 1 / 2 /R, for a sample 
of HCs and for 123385. The luminosity L of the central 
star is given by the IRAS luminosity for HCs without an 
embedded Hn regions, and it represents an upper limit. 
For those with an embedded Hn region, the luminosity of 
the central star has been deduced from the radio contin- 
uum flux of the embedded Hn region (Wood & Churchwell 

HSS9- 

The drawn line indicates the steepest power law, T tx 
(L 1 / 2 /R) - 5 . It is evident that it is not possible to match 
our temperature estimate with that expected for such a 
power law. To solve this problem we examine three possi- 
bilities: 

1) We might be overestimating the source luminos- 
ity. In Sect. 14.11 we discussed that the luminosity of the 
embedded source is of the order of 10 4 solar luminosi- 
ties. This was based on the hypothesis that the bulk of 
the IRAS fluxes at 60 and 100 fjm is due to 123385. 
Unfortunately, the IRAS angular resolution is not suffi- 
cient to confirm this hypothesis. Hence, the "bright ring" 
surrounding 123385 might emit a consistent fraction of 
the flux at 60 and 100 /xm. In that case, the luminos- 
ity of 123385 might be lower. A temperature of 42 K 
at 0.017 pc in an optically thick core implies a central 
source of ~ 1.5 1O 3 L . Altough well below the value pre- 
dicted in Sect. l3~3~2l this luminosity is consistent with an 



intermediate-mass protostar deriving its luminosity from 
accretion, as shown by Behrend & Maeder (2001 ). In fact, 
the authors performed calculations of the evolution of pro- 
tostars with mass from 1 to 85 Mq , assuming growing ac- 
cretion rates. For an accretion luminosity of ~ 1.5 10 3 Lq, 
the protostar mass is ~ 6 Mq. 

2) Another possibility is that we have underestimated 
the core diameter. This is possible if the source is not 
gaussian: in fact, if the core is spherical, the diameter es- 
timated in Table requires a correction (see Panagia & 
Walmslev ll978|l . In our case one obtains ~ 2"(0.047 pc), 
for the diameter of the core. This brings the luminosity 
to 3 10 3 L Q , still much less than the value obtained in 
Sect. |4~T1 from the continuum spectrum. 

3) Finally, it is possible that the temperature profile, 
is steeper in case of a non-spherical geometry, as in the 
presence of a disk. Cesaroni et al. I|1998J) have found that 
HCs are likely to be disk-like structures, in which the tem- 
perature profile is T oc R~°' 7B . With this power-law it is 
possible to obtain a lower temperature at 0.017 pc with a 
luminosity of 1.6 10 4 L Q . 

Which solution is correct? Concerning the second hy- 
pothesis (underestimate of the core radius), we have al- 
ready seen that it is very unlikely. This is also valid for 
the third one (disk-like structures): in fact, to measure 
such a low temperature, the disk should lie along the line 
of sight, but unfortunately this is more likely "face-on", 
because the outflow detected by Molinari et al. ( 1998.EJ) lie 
approximately along the line of sight. The first hypothesis 
(overestimate of the core luminosity) is the most likely, 
but presently it is impossible to give a reliable answer, be- 
cause the correct source luminosity can be derived only by 
means of a map at ~ 100 /im at high angular resolution, 
which is not available. However, as already discussed, a 
luminosity as low as 10 3 Lq is consistent with a protostar 
in the accretion phase, whose mass, at present, is ~ 6 Mq. 



6. Conclusions 

We have used the IRAM-30m telescope, the Plateau de 
Bure Interferometer and the Very Large Array to observe 
the massive protostar candidate IRAS 23385+6053 in the 
C 18 (1-0), C 17 (1-0) and (2-1), CH 3 C 2 H (6-5), (8- 
7) and (13-12), NH 3 (1,1) and (2,2) lines. The following 
results have been obtained: 

— The source consists of two regions: a compact molec- 
ular core with a diameter ~ 0.03 + 0.04 pc and a halo 
which has a diameter of up to ~ 0.4 pc; the core has 
a kinetic temperature of Tk ~ 42 K and a H 2 volume 
density of the order of 7Jh 2 ~ 10 7 cm" 3 ; the halo has 
an average temperature Tk ~ 15 K and an average H 2 
volume density uh 2 ~ 3 10 5 cm" 3 ; 

— based on the continuum spectrum, the source lumi- 
nosity can vary between 1.6 10 4 and 1.5 10 2 Lq: we 
believe that the upper limit is closer to the real source 
luminosity; 
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Fig. 20. Kinetic temperature measured at the surface of 
molecular cores plotted against the distance-independent 
ratio L 1 / 2 /9, where 8 is the angular diameter of the core, 
and L the luminosity of the embedded (proto)star. Open 
circles indicate HCs (Kurtz et al. 2000) while the filled cir- 
cle corresponds to 123385. The straight line represents the 
steepest power law for spherical cores, T oc (X 1 / 2 /9d) - 5 , 
expected for a dusty core heated by an embedded star 
(Doty & Leung RMQ . 



— by comparing the H2 volume density derived at differ- 
ent radii from the IRAS source, we deduce that the 
halo has a density profile of the type nn 2 oc r~ 23 , sug- 
gesting that it might be gravitationally unstable. This 
is further supported by the fact that, assuming a den- 
sity profile as above, the virial mass is much lower than 
the gas mass. 

— On the basis of the steepest theoretical temperature 
profile expected in cores heated by embedded stars 
(T oc i? -0 5 ), it is impossible to explain a temperature 
of 42 K at the core radius, if the source luminosity is 
L ~ 10 4 L Q . We propose three possible solutions: 1) 
we are overestimating the stellar luminosity. A tem- 
perature of 42 K at the core radius (~ 0.017 pc) is 
possible if the stellar luminosity is ~ 10 3 L Q . This solu- 
tion is consistent with an intermediate-mass protostar 
deriving its luminosity from accretion; 2) we are un- 
derestimating the core radius in case of non-gaussian 
source, but this solution is unlikely because the correc- 
tion that we should apply is very small; 3) the source 
is not spherical: in that case the temperature profile is 
expected to be steeper (T oc R~ - 75 ) and it is possi- 
ble to reach a temperature as low as 42 K at the core 
radius even with luminosities as high as 10 4 Lq. 

The first hypothesis seems the most likely, although 
the real source luminosity can be derived only with 
high resolution maps at FIR wavelengths, which will 
become feasible only with the HERSCHEL satellite. 
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Fig. 9. Left panel: NH3(1,1) integrated map over the velocity range (—48, —52) km s -1 . Contour levels range from 
0.012 (~ 3<t) to 0.024 by 0.004 Jy beam -1 . The channel spacing of the integrated spectrum (shown over the map) 
is ~ 1.13 km s -1 . The grey-scale represents the C 18 (1-0) main line emission (Fig. EJ|. The cross indicates the core 
position detected with the PdBI. Right panel: same as Left for NHs(2,2). Contour levels range from 0.012 to 0.025 by 
0.004 Jy beam -1 . The channel spacing is ~ 4.9 km s . The ellipses in the bottom-right represent the VLA synthesised 
beam. 
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Fig. 19. Left panel: map of the C 17 (1-0) line in the peak channel (at —50.5 km s -1 ). The solid contour represents 
the observed 2a level (~ 0.024 Jy beam -1 ), while the dashed contours indicate the emission of a "core-halo" source 
model, as described in Fig. El Levels range from 0.024 to 0.14 by 0.024 Jy beam -1 . Right panel: same as left for a 
source model without the central core. Both contours represent the 2a observed level. 



